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INTRODUCTION 


This  report  supplies  a  program,  with  examples,  for  the  analysis  of  rectangular  microstrip 
antennas.  The  formulas  upon  which  the  program  is  based  are  also  provided.  The  theory  from 
which  these  formulas  were  obtained  is  based  on  the  "cavity  model"  of  the  microstrip  antenna 
developed  at  the  University  of  Illinois  by  Lo,  Richards,  et  al.  Details  of  the  theory  can  be 
found  in  the  references  listed  in  the  bibliography  at  the  end  of  the  report 

The  first  chapter  contains  the  formulas  and  the  definitions  of  the  electrical  and  geometri¬ 
cal  parameters  used  in  the  program.  Chapter  two  lists  the  FORTRAN  program  which  imple¬ 
ments  these  formulas.  Chapter  three  contains  examples  of  the  program’s  use.  It  includes  an 
example  which  illustrates  the  use  of  the  two-port  analysis  feature  of  the  program  to  determine 
the  "tuning  range”  that  a  variable  capacitor  loading  one  port  would  have  to  have  in  order  for  the 
radiator  to  produce  any  polarization  from  left  hand  circular  to  right  hand  circular  polarization. 


3. 


CHAPTER  1:  FORMULAS  AND  DEFINITIONS 


For  completeness,  the  basic  formulas  and  procedures  used  to  compute  the  pattern  and 
impedance  of  a  microstrip  antenna  are  reviewed.  The  geometrical  parameters  used  in  the  for¬ 
mulas  that  follow  are  defined  in  figure  1. 


I  GENERAL  NOTATIONS 


The  following  notations  are  used  in  the  formulas  within  this  appendix. 

(1)  k0  —  2ir f/c  where  /is  the  frequency  and  c  is  the  speed  of  light  in  free  space. 


(2)  t,0-  377O. 


(3)  «0m  —  1  for  m  —  0  and  2  otherwise. 

(4)  (r,9,<t>)  is  the  coordinate  of  a  point  in  spherical  coordinates.  The  direction  perpendicu¬ 
lar  to  the  ground  plane  corresponds  to  9  -  0.  The  line  9  -  ir/2,  —  0  is  the  x  axis 

while  4>  —  ir/2  is  the  y  axis. 


(5)  pm  -  (k2  -  (mn/a)2)'*  (the  branch  is  irrelevant). 


(6)  k  -  where  the  branch  is  also  irrelevant,  «,  is  the  relative  dielectric  con¬ 

stant  of  the  dielectric  substrate,  and  6  is  the  loss  tangent  of  the  dielectric  substrate 
(later  to  be  replaced  by  the  "effective  loss  tangent."  See  V.) 


(7)  *«> 


j'/i 

■^5L  cosIpm(f»-^)lcos(mirx/a),  for  y  >  yu  and 
a 


*0m 

a 


'h 

cos(p„y)cos(mirx/a),  fory  <  yx. 


(8)  jQ(x)  -  sin(x)/x  (the  spherical  Bessel  function  of  zero  order). 

(9)  A  is  the  skin  depth. 


II  PATTERN  AND  RADIATED  POWER 

Radiated  power,  Pm ,,  is  computed  in  subroutine  VRAD.  This  routine  calls  the  double 
integration  routine,  VDOUBL,  which  applies  4-point  Gaussian  quadrature  recursively  to 
integrate  the  power  pattern  supplied  by  VPPAT.  The  function  VPPAT  calls  VP  AT  which  com¬ 
putes  the  complex  polar  pattern  of  the  antenna  by  application  of  the  following  formulas: 

F  Jkptyob  "  g0mcos(rnffX|  /a)  |  afco*  A 

r  2m  pmbnn(pmb)  °(  2 a  Jl  * 

■  |i[co.W, -  coslA,(»-j.|)l] 


2 


Figure  1.  Geometry  and  Idealized  Feeds  for  the  Rectangular 
Microstrip  Antenna 
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-  yQpmbs\n(pmb)eJko>l's'T0s'T*‘  +  yAro6sintfsin0 

•  [cosW,)/0^  -  cosl/»m(*— >,))]] 

•  |(pmb)2  -  (fco*sin0sin^)2J  | 

Ee  —  ko(-  Fxsin4>  +  Fycos<t>) 

E*  “  — ko(Fxcostj>  +  Fysin(>)cos0 

where  Ee  and  E+  are  the  9  and  <f>  components  of  the  electric  far  field,  Fx  and  Fy  are  the  jrand  y 
components  of  F,  and  Jr  and  y  are  the  unit  vectors  in  the  x  and  y  directions. 


III  STORED  ENERGY  AND  OHMIC  LOSSES 


The  stored  electric  energy,  is  computed  from 


4irfWE  -  er*0f  <M)^o  £  1 - 1 - — 

«-o  |pmAsin(pm6) 


r 

y\  <I>*,(x1^1)|V(pB,yi)  +  (b-.y,)|<I><2)(x1,.y,)|  N[pm(b-yi)]  , 


where 

N(z)  -  y[/00’21mz)  +  y0(2Rez)l. 

The  dielectric  loss  is  found  from 
P-4nfbWE 

The  copper  loss  is  determined  using 

^cu  __  A 
Pi  bt 

(at  resonance).  All  these  quantities  are  computed  within  subroutine  VENLS. 

IV  IMPEDANCE 

The  impedances  are  computed  in  subroutines  VZ1  and  VZ2.  The  former  is  called  by  VZ2 
to  compute  zti  and  Z22  while  Z12  is  computed  within  VZ2  by  the  following  formula  for  >2  > 


zi2"  -JkotVoT  ^-cosimnxi  / a)cos(mirxj / 

a  I  2o 


4 


cosMft-^lcosO^i)  1 
pmsin(pmb)  (• 

For  >1  —  y2,  this  series  is  accelerated  by  writing  it  as 
|  cos(A>'i)cosU(6-.y1)l 


Z12  “  -jk0tT)0\ 


■ 


kasin(kb) 

+  ^,^-cos(mirxi/a)cos{nmx2/a)ji\-^^- 

£ \a  (  2fl 

cos(/>m>>,)cos[pm(ft-yi)] 


pmsin(p„b) 


ar 

mv 


jkotyc? 


7T 


n(x t+x2)  I  ,  _(  ir(x,-x2) 


1 


+  F 


ir(x\+x2+d) 


+  F 


ir(x]+x2-d) 


+  F 


it  (x\—x2+d)  I  |  ir(x\-x2-d) 

a  J+  |  a 


where  r  —  1  for  b  >  j>!  >  0  and  r  —  2  for  —  0  or  >]  —  b.  The  driving  point  impedance,  zt), 
is  computed  using  the  accelerated  formula  for  zt2  with  x2  and  y2  replaced  by  xi  and  >>i,  respec¬ 
tively.  Similarly,  z12  is  computed  with  x(  and  y\  replaced  by  x2  and  y2. 

The  function  F(x)  is  related  to  Clausen’s  integral  and  is  given  by 

(mx) 


F(x)  -  Xcos- 


m- 1 


m 


This  function  is  written  in  terms  of  Inx  and  a  rapidly  converging  series  of  Chebyshev  polynomi¬ 
als.  It  is  evaluated  in  function  VF. 


V  EFFECTIVE  LOSS  TANGENT 

The  "effective  "  loss  tangent  is  found  by  first  computing  the  fields  within  the  "cavity" 
based  on  k  found  from  the  actual  loss  tangent  of  the  substrate.  From  these  fields,  computa¬ 
tions  of  the  electric  stored  energy,  the  radiated  power,  and  the  power  loss  in  the  dielectric  and 
copper  are  made.  From  these  quantities,  the  antenna "  Q "  is  computed  from 

E 

PnA  +  F<i  +  f*Cu 

An  "effective"  loss  tangent,  6eff,  is  defined  as 


This  loss  tangent  is  then  used  to  compute  an  improved  k  and  the  whole  process  is  repeated  to 
find  new  (and  more  accurate)  predictions  of  the  stored  energy  and  losses.  A  new  fie(r  is  found, 
and  so  on.  The  program  as  supplied  computes  a  twice  iterated  6eff.  However,  for  thin  sub¬ 
strates,  the  procedure  converges  after  a  single  iteration  and  the  first  8eff  computed  is  adequate. 
A  simple  modification  of  the  program  will  eliminate  the  second  iteration. 


r 

^>**v 


CHAPTER  2:  PROGRAM  LISTING 


Hie  FORTRAN  program  listed  in  this  chapter  was  implemented  on  the  CYBER  17S  com¬ 
puter  located  at  the  University  of  Illinois,  Urbana,  IL.  The  program  uses  CDC’s  "extended 
FORTRAN"  and  the  Graphics  Compatibility  System  (GCS)  produced  by  the  United  States  Mili¬ 
tary  Academy.  Names  of  subroutines  provided  by  GCS  all  begin  with  the  letter  "U."  These 
GCS  routines  are  used  in  certain  input/output  subroutines  including  those  that  plot  results. 
Such  routines  have  been  listed  below  in  a  section  entitled  "INPUT/OUTPUT  AND  PLOT¬ 
TING."  It  is  this  section  of  die  program  which  is  rather  strongly  installation  dependent  and 
would  probably  require  user  modification. 

The  other  major  sections  of  the  program  listed  below  are  the  "MAIN  PROGRAM"  and 
the  "NUMERICAL"  sections.  The  former  controls  the  flow  of  execution  of  the  program  while 
the  latter  computes  the  impedance,  pattern,  etc.  Both  these  sections  are  fairly  transportable, 
particularly  the  NUMERICAL  section.  Only  a  few  non-ANSI  FORTRAN  statements  and  rou¬ 
tines  are  used  in  these  sections  and  these  can  be  easily  eliminated  or  modified. 

The  overall  simplified  flow-chart  of  the  program  is  shown  in  fig.  2.  Other  details  of 
specific  subroutines  can  be  obtained  by  referring  to  their  respective  documentation  in  comment 
cards. 
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Input  Parameter* 

[viNi]vs*trr*o 


Compute  P<j,  Pcu 

Wg 

[vENLS] 


Compute  Prod 
[VRAD] 


Compute  Q  ,  8«ff 
[vmain] 


IT.IT-H 

[vmain] 


Compute  81  Plot  Z,n 
[VZ1]  [VPLTZ] 


Port2\  Yes 
Loaded  / — - — 

v  ?  / 


Compute  8  Plot 
s  Parameter* 
[VZ2][VPLTS] 


Change  Frequency 
[UMAIN] 


Compute  8  Plot 
Zin 

[vze]  [vpltz] 


Find  Current  in 
Port  2 

[VMAIN]  [VLOAD] 


Plot  Pattern 
[VPLTPT] 


Figure  2.  Simplified  Flow  Chart  for  the  Program 
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MAIN  PROGRAM 


PROGRAM  VMAIN( INPUT, OUTPUT .TAPE1=0UTPUT,TAPE2=INPUT, RESULT,  N/S 

1  TAPE3=RESULT) 


This  is  the  main  program  of  a  group  of  routines  that  computes 

(1)  Input  impedance  of  a 

(a)  Single  port  rectangular  microstrip  antenna: 

(b)  Two  port  rectangular  microstrip  antenna  with  one  of  its 
ports  loaded  by  a  specified  impedance  (in  subroutine 
"VLOAD" ) ; 

(2)  The  "s"  parameters  of  a  two  port  microstrip  antenna; 

(3)  The  radiation  pattern  of  a 

(a)  Single  port  rectangular  microstrip  antenna: 

(b)  Two  port  rectangular  microstrip  antenna  with  one  of  its 
ports  loaded  by  a  specified  impedance. 


REFERENCES:  The  method  used  is  described  in  the  following 
publications: 

[1]  Y.  T.  Lo,  D.  Solomon.  W.  F.  Richards,  "Theory  and  Experiment  on 
Microstrip  Antennas,*  IEEE  TRANS.  ANTENNAS  PROPAGAT.  Vol. 
AP-27,  pp.  137-145,  MAR  79. 

[2]  Y.  T.  Lo,  W.  F.  Richards,  D.  D.  Harrison,  "An  Improved  Theory 
for  Microstrip  Antennas  and  Applications,"  RADC-TR  INTERIM 
REPORT  (PART  I),  DEC  78. 

[3]  W.  F.  Richards,  Y.  T.  Lo,  D.  D»  Harrison,  "Improved  Theory 
for  Microstrip  Antennas,*  IEE  ELECTRONICS  LETTERS ,  Vol.  15, 
pp.  42-44,  JAN  79. 

LIMITATIONS: 

The  current  version  does  not  include  an  estimate  of  surface  wave 
power  as  this  computation  is  currently  under  critical  evaluation. 
Thi3  version  also  requires  the  specification  of  the  so  called 
"effective  feed  width."  This  parameter  arises  from  an  attempt  to 
idealize  the  fields  in  the  vicinity  of  a  coaxial  or  microstrip  feed 
so  that  the  source  can  be  considered  as  a  uniform  current 
ribbon  of  width  D  (the  effective  width)  flowing  from  the  patch  to 
the  ground  plane.  Since  the  observed  shift  of  impedance  loci  into 
inductive  half  of  the  Smith  Chart  depends  rather  strongly  on  the 
field  distribution  in  the  vicinity  of  the  feed,  this  idealization 
needs  some  refinements  and  a  more  rigorous  treatment  of  this 
problem  is  under  way.  For  the  present,  the  user  should  try  some 
different  values  of  D  until  he  rinds  one  which  fits  measured  results 
most  closely.  The  representation  of  the  fields  for  frequencies  far 
away  from  resonance,  say  near  the  mean  of  two  widely  spaced  adjacent 
resonant  frequencies,  is  currently  not  sufficiently  accurate  for  all 
applications.  We  will  do  further  research  to  develop  better  compu¬ 
tations  in  this  regime. 

IMPLEMENTATION  REQUIREMENTS:  Except  for  the  input/output  which 
relies  heavily  upon  the  graphics  capabilities  of  the  GCS  system 
developed  by  the  United  States  Military  Academy,  the  program  is 
written  in  ANSI  FORTRAN  and  should  be  relatively  transportable. 

All  GCS  subroutine  names  begin  with  a  "U"  in  this  program. 

Some  of  the  input/output  utilizes  extended  features  or  CDC’s  FORTRAN 
as  implemented  on  the  University  of  Illinoi's  CYBER  175  (NOS  V.  4.7) 
and  will  have  to  be  modified  for  use  on  other  systems. 

Non-ANSI  FORTRAN  statements  are  flagged  as  N/S. 

USER  INSTRUCTIONS:  The  parameter  descriptions  and  options  are 
explaned  through  the  use  of  examples  provided  with  this  listing. 
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REAL  LOSS,  KO.  LOSSO 
INTEGER  P,  ANS,  PO 

COMPLEX  2,  211 f  212.  222,  2L,  2IN,  12 
COWON  /DELTA/  DELTA 
COWON  /OPT/  ANS 

COWON  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP ,  YP,  XPP,  YPP,  L, 

•  F6,  NFREQ,  KO,  TWOPI,  PD02A ,  PXPOA,  BOA,  BAND,  PI, 

•  ETAO.  K,  GAIN,  DELTA* 

COWON  /I2/  AO,BO,TO,DO,DlELO,LOSSO,SIGMO,PO,XPO,YPO,XPPO,YPPO, 

•  lo,foo,band6,deltao 

DATA  AO/O . / , B0/0 . / , TO /0 . / , DO /0 . / . DIELO/ 1 . / , LOSSO/O . / , SIGMO/580 . / , 

•  PO/1 / ,XP0/0 . / 'YPO/O . / ,XPPO/0 . / .YPPO/O . /,L0/"N"/,F0O/0 . / , 

•  BANDO^O./.DELfAO/O.O/, C/30000./, PI/3. 1415926535898/, 

•  TWOPI/6. 2831853071796/, ETAO/377./ 

CALL  USTART 

1  CALL  VIN1 
IT  a  0 

DELTA  =  LOSS 

Find  the  nearest  resonant  frequency  to  the  specified  center  frequency 


CALL  VSEARC  (M,  N) 

KO  a  SQRT(((M»h/A)«2  +  (N»PI/B)*»2)/DIEL) 

PD02A  *  PI  •  D  /  (2  •  A) 

PXPOA  =  PI  •  XP  /  A 
BOA  a  B  /  A 

Compute  "closed  form"  sum  of  asymptotic  expression  of  the  summand 
for  the  driving  point  impedance  series 

SI  a  VS(XP,  XP) 

IF  (P  .LT.  2)  GO  TO  2 
a  VS(y 


S2 


XPP,  XPP) 


S3  a  VS(XP,  XPP) 

Compute  stored 
loss,  PD,  for  " 
always  computed 


electric  energy,  WWE,  copper  loss,  PCU,  and  dielectric 
effective  loss  tangent",  DELTA.  (Dielectric  loss  is 
using  the  actual  loss  tangent,  LOSS). 


2  CALL  VENLS  (DELTA,  WWE,  PCU,  PD) 

Compute  radiated  power,  PRAD,  for  effective  loss  tangent,  DELTA 

CALL  VRAD  (DELTA,  PRAD) 

POWER  a  PRAD  +  PCU  +  PD 
Q  a  2  •  WWE  /  POWER 
DELTA  =  1/Q 
IT  a  IT  +  1 

Iterate  the  calculation  twice  to  ensure  proper  value  of  DELTA 
is  obtained. 

IF  (IT  .LT.  2)  GO  TO  2 

Compute  the  pattern  along  the  zenith  direction  to  determine 
antenna  gain. 

CALL  VPAT  (0.,  0..  211.  222.  DELTA,  (1.,  0.),  XP,  YP) 

E  a  211 *C0N JG (211 )+222»C0N JG ( 222 ) 

GAIN  a  10 .•ALOG(2"TWOPI*E/(ETAO*PRAD) ) 

1  /AL0G(10.) 

F  a  FO  -  (NFREQ/2)*DELTAF 
N02  a  NFREQ/2  +  1 
WRITE  (1,7) 

READ  (2,8)  ANS 
CALL  E0F(2) 

if  Tans  .eq.  o)  ans  a  3 

DO  3  K  a  1  NFREQ 
KO  a  TWOPI  •  F  /  C 


N/S 
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GO  TO  (10,  20),  P 

Find  the  driving  point  impedance  of  the  one  port 

10  CALL  VZ1  (Z,  DELTA,  XP,  YP,  SI) 

Input  data  to  the  plotting  program,  VPLTZ 

CALL  VPLTZ  (Z,  F,  0) 

GO  TO  30 

Compute  the  "z"  parameters  of  the  two  port 

20  CALL  VZ2  (Z11.  Z12,  Z22.  DELTA,  SI,  S2,  S3) 

IF  (L  .EQ.  "N")  GO  TO  202 

Compute  impedance  of  load  on  port  two  of  the  microstrip 
CALL  VL0AD  (F,  ZL) 

Compute  the  input  impedance  as  seen  at  port  one  of  the  loaded 
microstrip  antenna. 

ZIN  =  Z11  -  Z12«2  /  (Z22  +  ZL) 

CALL  VPLTZ  (ZIN,  F,  0) 


Compute  the  current  flowing  through  the  load  at  port  two. 
12  =  -Z12  /  (Z22  +  ZL) 


Depending  on  the  options  chosen,  compute  the  pattern  of  the  antenna. 

IF  (ANS  .EQ.  1  .OR.  (ANS  .EQ.  2  .AND.  K  .EQ.  N02)) 

1  CALL  VPLTPT  (12,  F) 

GO  TO  3 

Input  the  two  port  parameters  to  the  "s"  parameter  plotting 
program  for  the  ease  of  a  non-loaded  two  port  antenna. 


20 2  CALL  VPLTS  (F,  Z11.  Z12,  Z22.  0) 

30  IF  (ANS  .EQ.  1  .OR.  (ANS  .EQ.  2  .AND.  K 
1  CALL  VPLTPT  (12,  F) 

3  F  =  F  +  DELTAF 

IF  (P  .EQ.  1)  GO  TO  5 
IF  (L  .NE.  "N")  GO  TO  5 


.EQ.  N02)) 


Plot  "s"  parameters 

CALL  VPLTS  (F0,  Z11,  Z12,  Z22,  1) 
GO  TO  6 


Plot  the  input  impedance  to  the  microstrip. 

5  CALL  VPLTZ  (Z,  F0,  1) 

6  WRITE  (1.9) 

READ  (2,11)  ANS 
CALL  E0F(2) 

IF  (ANS  .EQ.  0)  ANS  =  "Y" 

IF  (ANS  .EQ.  "Y")  GO  TO  1 
STOP 

7  FORMAT  ("Choose  an  option:",/, 

1  "  (1)  Plot  patterns  at", 

2"  all  frequencies.",/, 

3  "  (2)  Plot  pattern  only", 

4"  at  center  frequency.",/, 

"  (3)  Plot  no  patterns",/, 

tion  (1 ,  2,  or  3):  ",) 


9  FORMAT  ("Continue  (type  Y  or  N)  ",) 


N/S 
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NUMERICAL 


SUBROUTINE  VDOUBL  (F,  M,  A,  B,  C,  D,  INT) 


PURPOSE:  This  subroutine  performs  the  integral  from  0  to  A  of  the 
integral  from  0  to  B  of  the  externally  declared  function 
F(X ,  i) . 

PARAMETERS:  The  parameters  are  as  stated  above  and  M  is  the  log  base 
two  plus  one  of  the  number  of  point  Gaussian  quadrature 
formula  used  in  the  mechanical  quadrature.  INT  is  the 
integral . 


REAL  INT 
DATA  W  / 


DATA  S  / 


,  W(31),  S(31) 

1 .0000(3000000000 

.36268378337836 

.10122853629038 

.16915651939500 

.09515851168249 

.09654008851473 

.09117387869576 

.07819389576707 

.05866409347854 

.03427386291302, 

.00701661000947/ 


.2ll6035507792. 
.75540440835500 
.98940093499165 
.23928736225214 
.50689990893223 


.9  HI  1151154527 
2»»(M-1 ) 

2*L1  -  1 
/2 


=  (B-A 


LI 
N1 
FI 
F2 

Fl»  =  (D+C)/2 
INT  =  0. 

DO  2  J  =  LI, 
SUM  =  0. 

DO  1  K  =  LI 


D-C3/2 

B+A)/2 


N1 


T1 
T2 
T 
T 
T1 
T< 

SUM 
INT  =  INT 
INT  =  INT 
RETURN 
END 


N1 


F? 


FI 
F2 
T1 
T2 
F3  -  T1 
F4  -  T2 
SUM+WCK 

-  ;  wU 
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.65214515486255 

.31370664587769 

.18945061045507 

.14959598881656 

.06225352393665 

.09563872007927 

.06765209300440 

.07234579410885 

.05099805926238 

.02539206530926 

.33998104358486 

.16343464249565 

.96026985649754 

.45891677765723 

.04830766568774 

.33186860228213 


.58771575724076 


.997263 


948/ 


.34785484513745 

.22238103445337 

.18260341504492 

.12462897125553 

.02715245941175 

.09384439908060 

.06331192422695 

.06582222277636 

.04283589802223 

.01627439473091 


94457502307323, 

14447196158280, 

42135127613064, 

66304428693022, 

84936761373257, 

96476225558751, 


FI 


)»^F(T3,T4)+F(T3,T6)+F(T5,T4)+F(T5,T6)) 
A  F2 
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SUBROUTINE  VENLS ( DELTA , ESE , DL , CL ) 


PURPOSE:  This  subroutine  computes  the  electric  stored  energy  and  the  1 
copper  and  dielectric  losses  at  resonance. 

PARAMETERS:  DELTA  is  the  effective  loss  tangent.  ESE  is  the  computed  ' 
electric  stored  energy  (at  the  nearest  resonance  to  FO  as 
determined  by  subroutine  VSEARC) .  DL  is  the  dielectric  loss  1 
computed  using  the  actual  loss  tangent.  CL  is  the  copper  loss4 
computed  as  a  proportion  of  the  dielectric  loss. 


REAL  MPI, MPI0A.MPI0A2.LT, JO, N.KO 
COMPLEX  K . SPB , P ,  PYP , PSPB ,PBMYP, CPYP , CPBMYP 

C0»<0N  /II/  A.  B,  T,  D,  EPS,  LT,  SIGM.IP,  XP,  IP,  XPP,  YPP,  L, 

•  FO,  NFREQ,  KO,  TWOPI,  PD02A ,  TPXPOA,  BOA,  BAND,  PI, 

•  ETA,  KK,  GAIN.  DELTAF 


)  *  KO 


u-in  •  ■  unxti  •  vuui 

K  =  CSQRT ( CMPLX ( EPS , -DELTA*EPS ) 

SD  =  SQRT(2/(1000«KO«ETA*SIGM)) 

BMYP=B-YP 
CPYP=CCOS(K*YP) 

CPBMYP = CCOS ( K * BM YP ) 

SPB=CSIN(K«B) 

FIsYP*  ( CABS ( CPBMYP ) ) * *2 
F2=K1YP«(  CABS  (CPYP ) )  ”2 

ESE=(F1*N(K*YP ) +F2*N (K*BMYP ) ) / ( CABS (K«SPB) ) «»2 
M1=0 

100  CONTINUE 
M1=M1+1 

MPI =M1«3. 141592654 
MPIOA=MPI/A 
MPI0A2=MPI0A**2 
P=CSQRT (K*K-MPI0A2 ) 

PSPB=P*CSIN (P*B) 

PYP=P«YP 

PBMYP=P*BMYP 

CPYP=CCOS(PYP) 

CPBMYP=CCOS (PBMYP) 

CF= ( JO ( Ml *PD02A)*C0S ( MPI »XP / A ) ) «»2 
F=2. /( CABS (PSPB) )**2 
F 1 =YP*( CABS ( CPBMYP ))**2 
F2=BMYP*( CABS( CPYP ) ) **Z 
T3=(F1 »N (PYP )+F2»N (PBMYP ) ) *F 
IF  (M0D(M1 ,2)  .EQ.  0)  GO  TO  150 
SUBTOT  =  CF*T3 
GO  TO  100 

150  SUBTOT  =  SUBTOT  +  CF*T3 
ESEsESE+SUBTOT 

IF  (SUBTOT/ESE  .LT.  0.0001)  GO  TO  200 
GO  TO  100 

200  ESE=ESE»EPS*T*K0«3#ETA/(2«A) 

DL=2*LT*ESE 

CL=SD»DL/(T*LT) 

RETURN 

END 
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FUNCTION  VF(Z) 


PURPOSE:  This  routine  computes  a  sum  related 

to  the  integral  of  Clausen's  integral: 

Stan  from  k  =  1  to  infinity  of  cos(kZ)/k**3. 

METHOD:  A  Tchebyshev  series  was  derived  from 

expansions  given  in  Abromowitz  &  Stegun 
ana  is  summed  by  Clenshaw's  algorithm. 


REAL  Cl (9),  C2(9) ,  LN202 
DATA  Cl/  .06905464409766968, 

1  .00052071642034163, 

.00000505847024742, 

.00000006775156996, 

.00000000107462004, 

.00000000001896987, 

.00000000000036087, 

.00000000000000726, 

.00000000000000015/ 

DATA  C2  /. 2085520596900512 , 

.0029753216441284, 

.0000406999817326, 

.0000006926756791, 

.0000000132725892, 

.0000000002743922, 

.0000000000059818, 

.0000000000001356, 

.0000000000000032/ 

TWOPI  =  6.283185307179586 
PIBY3  =  1.047 19755 119b598 
ZETA3  =  1.202056903159594 
LN202  =  .3465735902799726 
X  =  ABS(Z) 

IX  =  X  /  TWOPI 
IF  (X  .LT.  0.)  IX  a  IX  -  1 
Y  =  X  -  TWOPI  *  IX 

IF  (Y  -GT.  3.141592653589793)  Y  =  TWOPI  -  Y 
KODE  a  1 

IF  (Y  .GE.  2.094395102393196)  KODE  =  2 
IF  (KODE  .EQ.  2)  Y  =  3.141592653589793  -  Y 
T  =  Y  /  PIBY3 

IF  (KODE  .EQ.  1)  T  =  T  /  2. 

T2TSM1  =  2.  •  (2.  •  T”2  -  1.) 

G1  =  0. 

G2  =  0. 

DO  4  L  a  1,  9 
FACTOR  a  T2TSM1 

IF  (L  .EQ.  9)  FACTOR  a  .5  •  FACTOR 
GO  a  FACTOR  •  G1  - 
GO  TO  (1.2).  KODE 

1  C  a  CI(IO-L) 

GO  TO  3 

2  C  a  C2(10-L) 

3  GO  a  GO  +  C 
G2  a  G1 
G1  a  GO 

4  CONTINUE 
GO  a  GO  •  (Y 
GO  TO  (5.6), 

5  VF  a  ZETA3 

IF  (Y  .NE.  0.)  VF  a  VF  +  Y”2  •  (0.5  •  ALOG(Y)  -  0.75  -  '• 
GO  TO  7 

6  VF  a  -0.75  •  ZETA3  ♦  Y»*2  •  (LN202  -  2.  •  GO) 

7  RETURN 
END 


G2 


/  TWOPI )»»2 
KODE 


.»G0) 
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* 


SUBROUTINE  VL0AD(F,  ZL) 


PURPOSE: 


This  is  an  example  of  the  format  of  user  supplied  subroutine 
VLOAD.  The  purpose  of  the  subroutine  is  to  return  a  load 
impedance,  stored  in  ZL,  at  a  frequency,  F.  Many  of  the 
electrical  and  geometrical  parameters  or  the  antenna  are  ' 


electrical  and  geometrical  parameters  of  the  antenna  are 
available  to  this  routine  through  common  block  /II/. 

Also,  the  parameter  nKn  is  availble  through  this  common 
block.  This  parameter  is  a  "DO"  index  for  the  loop  that 
increments  frequency  in  the  main  program.  Thus,  as  was  done 
in  this  example,  and  "IF"  statement  testing  to  see  if  K  is 
one  or  not  can  be  included  so  that  data  can  be  input  by  this 
program.  Another  common  block,  /LDI D/  contains  a  string 
which  allows  one  to  give  a  verbal  description  of  the  type  of 
load  defined  in  the  subroutine.  The  description  must  he  40 
characters  or  less  in  the  CYBER  system. 


COMPLEX  ZL 
INTEGER  STRING(4) 

C0JW0N  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP,  YPP,  L, 

•  FO,  NFREQ,  KO,  TWOPI,  PD02A,  PXPOA,  BOA,  BAND,  PI, 

•  ETAO,  K.  GAIN,  DELTA* 

COMMON  /LDID/  STRING 

IF  (K  .GT,  1)  GO  TO  1 

PRINT  *,  "INPUT  CAPACITANCE  IN  PICOFARADS:  C  =  ",  N/S 

READ  •,  C 

ENCODE  (40,  2,  STRING)  C  N/S 

c  -  c  *  1  E-6 

1  ZL_=  CMPLX(0.,-1./(F»TW0PI*O) 

2  FORMAT  ("CAPACITIVE  LOAD:  C  a  ",E8.2,"  ") 

RETURN 

END 
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PURPOSE:  This  routine  returns  a  load  impedance  of  zero,  (the 

impedance  of  a  short).  It  is  loaded  as  the  defualt  VLOAD 
subroutine. 


SUBROUTINE  VLOAD(F,  ZL) 

INTEGER  STRING (4) 

COMPLEX  ZL 

COMMON  /LDID/  STRING 

DATA  STRING  /"short  circuit  "/ 

ZL  =  (0.  ,0. ) 

RETURN 

END 
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PURPOSE:  This  group  of  functions  compute  certain  quantities  used  in  the1 
evaluation  of  electric  stored  energy  computed  in  subroutine  1 
VENLS.  .  .  1 


REAL  FUNCTION  M  (Z) 

COMPLEX  Z 
REAL  JO,  10 

M  =  (I0(2»AIMAG(Z))  -  J0(2»REAL(Z)) )/2 

RETURN 

END 

REAL  FUNCTION  N  (Z) 

COMPLEX  Z 
REAL  JO.  10 

N  =  (J0(2*REAL(Z) )  +  I0(2»AIMAG(Z)))/2 

RETURN 

END 

REAL  FUNCTION  10  (X) 

T  =  EXP(X) 

10  =  1. 

IF  (X  .NE.  0.)  10  =  (T  -  1/T)  /  (2*X) 

RETURN 

END 


nno 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


SUBROUTINE  VPAT  (THETA,  PHI,  ETHETA,  EPHI,  DELTA,  I,  X,  Y) 


*  PURPOSE:  To  evaluate  the  electric  far  field  at  direction 

»  (THETA, PHI) 

*  PARAMETERS:  ETHETA  and  EPHI  are  the  THETA  and  PHI 

*  components  of  electric  field.  DELTA  is  the 

*  loss  tangent  used  in  the  computation. 


COMPLEX  ETHETA,  EPHI,  FX,  FY,  PM,  KSQ,  PMB,  PMY,  CPY,  CPBMY. I, XF. 

*  F,  PMBSPB,  EKYSS,  EKBSS,  XT,  YT,  XTERM ,  YTERM,  D2,  FACTOR 
REAL  KO,  MPI ,  KASC,  KBSS,  KYSS,  JO 

COMMON  /II/  A,  B,  T.  D,  DIEL,  LOSS.  SIGM,  P,  XP,  YP,  XPP ,  YPP,  L, 

*  FO,  NFREQ.  KO,  TWOPI,  PD02A ,  TPXPOA,  BOA,  BAND,  PI, 

*  ETAO,  K,  GAIN,  DELTAF 
PXOA  a  PI  »  X  /  A 

FX  =  (0.  ,0.) 

FY  =  (0.,0.) 

M  =  0 

CT  =  COS (THETA) 

ST  =  SIN (THETA) 

CP  =  COS (PHI) 

SP  =  SIN (PHI) 

KO  •  Y  *  SP  •  ST 

KO  *  B  *  SP  *  ST 

KO  *  A  *  CP  •  ST 


KYSS 

KBSS 

KASC 


KSQ  =  DIEL  *  K0**2  *  CMPLX(1.,  -DELTA) 
1  MPI  a  M  •  PI 

PM  =  CSQRT (KSQ  -  (MPI/A)**2) 

PMY  =  PM  *  Y 
PMB  =  PM  *  B 
CPY  a  CCOS(PMY) 

CPBMY  =  CCOS ( PMB-PMY ) 

PMBSPB  =  CSIN(PMB)  •  PMB 
EKYSS  =  CEXP(CMPLX(0 . , KYSS) ) 

EKBSS  =  CEXP(CMPLX(0.,KBSS)) 

D1  =  MPI**2  -  KASC** 2 
D2  =  PMB**2  -  KBSS**2 
XF  =  CPY  *  EKBSS  -  CPBMY 
F  =  CEXP ( CMPLX ( 0 . . KASC ) ) 

IF  (M  .NE.  2*(M/2))  F  =  -F 
F  =  F  -  1. 

IF  (ABS(DI)  .GT.  1.E-5)  GO  TO  2 


»>- 


->  Find  limiting  value  of  expression  when  D1  - >  0. 


XT  =  XF 

IF  (M  .GT.  0)  XT  a  XT  /  2. 

GO  TO  3 

2  XT  =  XF  *  F  *  CMPLX(0 . .KASC)  /  D1 

3  YT  a  -BOA  »  F  *  (PMBSPB  •  EKYSS  +  CMPLX(0 . ,KBSS)*XF)  /  D2 
FACTOR  a  COS(M  *  PXOA)  *  J0(M  *  PD02A)  /  (PMBSPB) 

IF  (M  .GT.  0)  FACTOR  a  FACTOR  *  2. 

M  a  M  +  1 

IF  (M0D(M.2)  .EQ.  0)  GO  TO  30 
XTERM  a  FACTOR  *  XT 
YTERM  a  FACTOR  *  YT 
GO  TO  1 

30  XTERM  a  XTERM  +  FACTOR  •  XT 
YTERM  a  YTERM  +  FACTOR  *  YT 
FX  a  FX  +  XTERM 
FY  a  FY  +  YTERM 

IF  (SQRT (CABS (XTERM) **2  +  CABS(YTERM)**2)  .LT.  0.0001  * 

•  SQRT ( CABS (FX)**2  +  CABS(FY)**2))  GO  TO  4 
GO  TO  1 
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4  FACTOR  =  K0**2  *  T  *  ETAO  *  B  /  TWOPI 
ETHETA  =  (-FX*SP  +  FY*CP)  »  FACTOR  *  I 
EPHI  =  -(FX*CP  +  FY*SP)  »  CT  •  FACTOR  • 
RETURN 
END 


REAL  FUNCTION  JO(X) 

JO  =  1. 

IF  (X  .NE.  0.)  JO  =  SIN(X)  /  X 

RETURN 

END 
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FUNCTION  VPPAT  (THETA,  PHI) 


PURPOSE:  THIS  SUBROUTINE  COMPUTES  THE  POWER  PATTERN  (TIMES  THE  SIN  OF 
OF  THE  POLAR  ELEVATION  ANGLE,  THETA)  FROM  THE  COMPLEX  PATTERN 
COMPUTED  IN  SUBROUTINE  VPAT. 

PARAMETERS:  (THETA,  PHI)  IS  THE  DIRECTION  OF  OBSERVATION  IN  SPHERICAL 
COORDINATES. 


COMPLEX  ETHETA,  EPHI 
COMMON  /DELTA/  DELTA 

COMMON  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP.  YPP,  L, 

*  f6,  NFREQ.  KO,  TWOPI,  PD02A,  TPXPOA,  BOA,  BAND,  PI, 

«  ETAO,  K,  GAIN,  DELTAF 

CALL  VPAT  (THETA,  PHI.  ETHETA.  EPHI,  DELTA,  (1.,0.)f  XP.  YP) 

VPPAT  =  SIN (THETA)  »  (ETHETA  *  CON JG( ETHETA)  +  EPHI  •  CONJG(EPHI)) 

RETURN 

END 
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SUBROUTINE  VRAD  (DELTA,  PRAD) 


PURPOSE:  To  evaluate  the  power  radiated  by  the  microstrip  antenna 

PARAMETERS:  DELTA  is  the  effective  loss  tangent,  and  PRAD  is  the 
radiated  power  computed  using  numerical  quadrature. 


COMMON  /II/  A,  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP.  XPP,  YPP,  L, 

•  FO,  NFREQ,  KO,  TWOPI,  PD02A,  TPXPOA,  BOA,  BAND,  PI, 

•  ETAO ,  K,  GAIN,  DELTA? 

EXTERNAL  VPPAT 

CALL  VDOUBL  (VPPAT,  3,  0.,  PI/2.,  0.,  TWOPI,  PRAD) 

PRAD  =  PRAD  /  ETAO 

RETURN 

END 
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FUNCTION  VS  (XI,  X2) 


PURPOSE:  THIS  FUNCTION  EVALUATES  THE  CONTRIBUTION  DUE  TO  THE  FIRST 

TERM  IN  THE  ASYMPTOTIC  SERIES  OF  THE  SUMMAND  IN  THE  Z-PARAM 
EXPRESSIONS.  (THIS  IS  USED  TO  APPLY  RUMMER'S  TRANSFORMATION 
TO  ACCELERATE  THE  CONVERGENCE  OF  THE  SERIES). 

PARAMETERS:  XI  AND  X2  ARE  THE  ABSCISSA  OF  THE  LOCATIONS  OF  PORTS  1 
AND  2,  RESPECTIVELY. 


COMMON  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP,  YPP,  L, 
*  FO.  NFREQ.  KO,  TWOPI,  PD02A ,  TPXPOA ,  BOA,  BAND,  PI, 

»  ETAO,  K,  GAIN,  DELTAF 

n-r  M  /VJ  iVn  \  1  I  •  r 


FI  =  PI 
F2  =  PI 
F3  =  PI 
F4  =  PI 
F5  =  PI 
F6  =  PI 
VS  =  VF 
VS  =  VF 
VS  =  - 
RETURN 
END 


(XI  +  X2  +  D) 


»  (XI  +  : 

•  (XI  -  : 

•  (xi  -  : 

(F3)  +  VF 
(FI)  +  VF 
(A/D)**2 


X2  -  D 
X2  +  D 
X2  -  D  J 


(F4)  +  VF(F5)  +  VF(F6) 
(F2)  -  0.5  *  VS 
•  (1/PI)**3  »  VS 


ooooonanoo 


SUBROUTINE  VSEARC  (MO,  NO) 


T 


PURPOSE:  This  subroutine  searches  for  the  combination  of  mode  indices, 
(MO, NO)  which  yields  the  resonant  wave  number  closest  two  the 
wave  number  of  free  space  at  the  chosen  center  frequency 
times  the  permitivity  of  the  dielectric. 


1 

2 

3 

4 

5 

6 

l 


REAL  MIN,  KG,  KMN 
INTEGER  V 

COJWON  /II/  A.  B.TT,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP,  YPP, 

•  f6.  NFREQ,  KO,  TSOPI.  PD02A,  TPXPOA,  BOA,  BAND,  PI 

•  ETAO,  K,  GAIN,  DELTA? 

DATA  C/30000./ 

KG  =  TWOPI  »  FO  •  SORT (DIEL )  /  C 
M  =  KG  •  A  /  PI 
N  =  0 

MIN  =  (M+1)  •  PI  /  A  -  KG 
MO  =  M  +  1 
NO  =  0 
V  =  1 

IF  (V  .EQ.  1)  GO  TO  2 
M  =  M  -  1 
GO  TO  3 
N  =  N  +  1 


KMN  =  SQRT ( ( M*PI /A ) **2  +  (N*PI/B)»*2) 
IF  (KMN  .LT.  KG)  GO  TO  4 

V  =  0 
GO  TO  5 

V  =  1 

T  =  ABS(KMN  -  KG) 

IF  (T  .GE.  MIN)  GO  TO  6 
MIN  =  T 
MO  =  M 
NO  =  N 

IF  (M  .NE.  0)  GO  TO  1 
IF  (V  .EQ.  1)  GO  TO  7 
N  =  N  -  1 
GO  TO  8 
N  =  N  +  1 

T  =  ABS(N*PI/B  -  KG) 

IF  (T  .GE.  MIN)  RETURN 

MO  =  0 

NO  =  N 

RETURN 

END 


L, 


23 


oooooooouooooo 


SUBROUTINE  V21  (Z,  DELTA,  X,  Y,  S) 


PURPOSE:  This  subroutine  computes  the  driving  point  impedance  of  a 
rectangular  microstrip  antenna  feed  at  point  (X,Y). 

PARAMETERS:  Z  is  the  complex  driving  point  impedance.  DELTA  is  the 
effective  loss  tangent,  (x, i)  is  the  coordinate  of  the  feed. 
S  is  the  "closed  form"  sum  of  the  asymptotic  form  of  the 
summand  for  Z.  (It  is  used  to  accelerate  the  convergence 
the  series. 


i 


COMPLEX  Z,  TERM.  K,  KSQ,  PM,  PMB,  PY,  PBMY ,  SUBTOT 
REAL  MPI,  KO,  LOSS,  JO 

COtfiON  /II/  A.  B.  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP.  YPP,  L, 

•  FO,  NFREQ,  KO,  TWOPI,  PD02A ,  PXPOA,  BOA,  BAND,  PI, 

•  ETAO,  KK,  GAIN,  DELTAF 
TAU  =  1.0 

IF  (Y  .EQ.  0.)  TAU  =  2. 

KSQ  =  DIEL  »  K0**2  *  CMPLX( 1 . , -DELTA) 

K  -  jpS^iPccOSCK^Y)  *  CC0S(K*(B-Y) )  /  (K«CSIN(K»B) ) 

M  =  0 

SUBTOT  =  (0.,  0.) 

PXOA  =  PI  •  X  /  A 

1  M  =  M  +  1 
MPI  -  M  *  PI 

PM  ='CSQRT(KSQ  -  (MPI /A) ”2) 

PMB  *  PM  •  B 
PY  =  PM  *  Y 
PBMY  =  PMB  -  PY 

TERM  =  (2*CC0S(PY)*CC0S (PBMY ) /(PM*CSIN (PMB) )  +  TAU»A/MPI)  • 

•  (COS(M»PXOA)  •  J0(M*PD02A))**2 
SUBTOT  =  SUBTOT  +  TERM 

IF  (3*(M/3)  .NE.  M)  GO  TO  1 
Z  =  Z  +  SUBTOT 

IF  (CABS (SUBTOT)  .LT.  0.001  *  CABS(Z))  GO  TO  2 
SUBTOT  =  (0.,0.) 

GO  TO  1 

2  Z  =  -Z  •  CMPLX(0. , 1 . )  •  KO  •  T  •  ETAO  /  A 
RETURN 

END 


oooooooooooooo 


SUBROUTINE  V22  (Z11,  Z12,  Z22 ,  DELTA,  SI,  S2,  S3) 


PURPOSE:  This  subroutine  computes  the  open  circuit  parameters  of  the 
two  port  with  port  one  at  (XP,YP)  and  port  two  at  (XPP, YPP) 
(where  these  parameters  are  in  common  block  II). 

PARAMETERS:  Z11.  Z12,  and  Z22  are  the  computed  open  circuit  parameters 
DELTA  is  the  effective  loss  tangent.  Si,  S2,  and  S3  are  the 
"closed  form"  sums  of  the  asymptotic  form  of  the  summands 
corresponding  to  Z11,  Z12,  and  Z22  summations,  respectively. 


COMPLEX  Z11.  Z12,  Z22.  TERM,  K,  KSQ,  PM,  PMB,  PYP,  PBMYPP , SUBTOT 
REAL  MPI,  KO,  LOSS,  JO 

COJWON  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P, XXP, YYP, XXPP, YYPP,  L, 

•  FO,  NFREQ,  KO,  TWOPI,  PD02A ,  PXPOA,  BOA,  BAND,  PI, 

•  ETAO,  KK.  GAIN.  DELTAF 
IF  (YYP  .LT.  YYPP)  GO  TO  10 

XP  =  XXPP 
YP  =  YYPP 
XPP  =  XXP 
YPP  =  YYP 
GO  TO  20 
10  XP  =  XXP 
YP  =  YYP 
XPP  =  XXPP 
YPP  =  YYPP 
20  TAU  =  1.0 

IF  (YP  .EQ.  0.)  TAU  =  2.  y 

KSQ  *  DIEL  •  K0«2  •  CMPLXO . , -DELTA) 

K  =  CSQRT(KSQ) 

KODE  =  1 

Z12  =  CCOS(K*(B-YPP))  *  CC0S(K*YP)  /  (K»CSIN(K»B) ) 

IF  (ABS(YP  -  YPP)  .LT.  0.001)  KODE  =  2 
IF  (KODE  .EQ.  2)  Z12  s  Z12  +  S3  *  A 
M  s  0 

SUBTOT  =  (0.,0.) 

PXPOA  =  PI  *  XP  /  A 
PXPPOA  =  PI  *  XPP  /  A 

1  M  =  M  +  1 
MPI  *  M  •  PI 

PM  =  CSQRT(KSQ  -  (MPI /A) ”2) 

PMB  =  PM  •  B 

PYP  =  PM  •  YP 

PBMYPP  =  PM  *  (B  -  YPP) 

TERM  =  2«CC0S( PBMYPP)  •  CCOS(PYP)  /  (PM*CSIN(PMB)) 

IF  (KODE  .EQ.  2)  TERM  =  TERM  +  TAU»A/MPI 

TERM  =  TERM  •  C0S(M»PXP0A)  •  COS(M»PXPPOA)  •  J0(M«PDO2A)”2 
SUBTOT  =  SUBTOT  +  TERM 
IF  (3*(M/3)  .NE.  M)  GO  TO  1 
Z12  =  Z12  +  SUBTOT 

IF  (CABS (SUBTOT)  .LT.  0.000001  *  CABS(Z12))  GO  TO  2 
SUBTOT  =  (0.,0.) 

GO  TO  1 

2  Z12  =  -Z12  *  CMPLX(0. , 1 . )  •  KO  •  T  •  ETAO  /  A 
CALL  VZ1  (Z11,  DELTA, XXP, YYP,  SI) 

CALL  VZ1  (Z22 ,  DELTA, XXPP, YYPP,  S2) 

RETURN 

END 
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SUBROUTINE  VZT0S(Z1 1 ,Z12 ,Z22 ,S1 1 ,S12 ,S22) 


PURPOSE:  This  program  converts  the  open  circuit  parameters  to 
scattering  parameters  refered  to  a  50  ohm  system. 

PARAMETERS:  The  open  circuit  parameters  are  Z11,  Z12,  and  Z22  and  are 
are  converted  to  the  scattering  parameters  $11,  S12,  and  S22. 


COMPLEX  Z11 ,Z12,Z22,S11 ,S12,S22 
Z0=50 

S11 =((Z1 1-Z0),(Z22+Z0)-Z12**2)/( (Z1 1+Z0)*(Z22+Z0)-Z12**2) 


S12=2,*Z0*Z12/( (Z1 1+Z0)*(Z22+Z0;-Z12**2)  „  M  v 

S22=(^Z11+ZO)»(Z22-ZO)-Z12««2)/((Z11+ZO)»(Z22+ZO)-Z12»»2) 


RETURN 

END 
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INPUT/OUTPUT  AND  PLOTTING 


SUBROUTINE  VIN1 


PURPOSE:  To  provide  input  of  antenna  parameters. 


INTEGER  P,  PO 
REAL  LOSS.  LOSSO 

COJMON  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP,  YPP,  L, 

*  FO.  NFREQ,  KO,  TWOPI,  PD02A,  TPXPOA,  BOA,  BAND,  PI, 

*  ETAO.  K.  GAIN,  DELTAF 

COMMON  /I2/  AO , BO , TO , DO , DI ELO , LOSSO , SIGMO , PO , XPO , YPO , XPPO , YPPO , 

*  lo \ fo6 , bAnd6 , deltAo 


WRITE  (1,1), 

1  FORMAT  (*  >  = 


i  a-a-  * 


«A#A#AIAiAIAiA$- 

•A»A#A#A»A# 

mAmAmAmAAAm 


A»A»n  / 

A»A»A4A4A»A«n  / 
A« A#~#A# A#A«n  ’/ ’ 
A*A*A*A*A*A*n  ’  /  ’ 
A»A#A«A#A»A«n(  ' 


OPT 


»- 


->  Erase  sequence  +  non-printing  buffer  characters 
for  Tektronix  graphics  terminals. 


WRITE  (1,2) 

FORMAT ( "  .  a  .  b  .  t  .  d  .diel .loss.sigm.p.") 

FORMAT  (F5.0fX.F5.0,X,F3.0,X,F3.0,X,F4.0,X,F4.0,X,F4.0,X,I1) 
READ  (2.3)  A,  B,  T,  D,  DIEL,  LOSS,  SIGM,  P 
IF(EOF(2) .EQ.O)  LOSS  =  LOSS/IOOO. 


T 

D 


A  =  AO 
B  =  BO 
=  TO 
=  DO 

DIEL  =  DIELO 
LOSS  =  LOSSO 
SIGM  = 

PO 

•  DIEL 


0. 

P  = 

•  D 


IF  (A  .EQ.  0. 

IF  (B  .EQ.  0. 

IF  (T  .EQ.  0. 

IF  (D  .EQ.  0. 

IF  (DIEL  .EQ.  0. 

IF  (LOSS  .EQ.  0. 

IF  (SIGM  .EQ. 

IF  (P  .EQ.  0) 

IF  (A  *  B  «  T 
AO  =  A 
BO  =  B 
TO  =  T 
DO  =  D 
DIELO  =  DIEL 
LOSSO  =  LOSS 
SIGMO  =  SIGM 
PO  =  P 

IF  (P  .GT.  1)  GO  TO  6 
WRITE  (1,4) 

FORMAT  (*  .  X’  .  y’  .") 
READ  (2.5)  IXP,  IYP 
FORMAT  (A5,X,A5) 

I  r  E0F(2) 

GO  TO  9 
WRITE  (1,7) 

FORMAT  (*  .  x’  .  y’ 

READ  (2.8)  IXP,  IYP 


OPT 

OPT 

OPT 

OPT 


OPT  N/S 


SIGMO 

*  LOSS  *  SIGM 


.EQ.  0.)  GO  TO  10 


.L.") 


nunv  |W/  xai  *  lit  *  IXPP,  IYPP,  L 

8  FORMAT  (A5,X,A5,X,A5,X,A5,X,A1) 

I  —  EOF (2) 

IF  (L  .NE.  "Y"  .AND.  L  .NE.  "N")  L  =  LO 
DECODE  (5, 12, IXPP)  XPP 
DECODE  (5, 12, IYPP)  YPP 
IF  (IXPP  .EQ.  "  ") 

IF  (IYPP  .EQ.  "  ") 

XPPO  =  XPP 
YPPO  =  YPP 
LO  s  L 

9  DECODE  (5, 12, IXP)  XP 


XPP  =  XPPO 
YPP  =  YPPO 


OPT  N/S 


OPT  N/S 

N/S 

N/S 


N/S 
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DECODE  (5,12 ,IYP)  YP 

IF  (IXP  .EC.  "  ")  XP  =  XPO 

IF  (IYP  .EQ.  "  ")  YP  =  YPO 

XPO  =  XP 

YPO  =  YP 

GO  TO  13 

10  WRITE  (1,11) 

11  FORMAT  INPUT  DATA  IS  BAD  /  JOB  ABORTED") 

STOP 

12  FORMAT  (F5.0) 

13  WRITE  (1,14) 

14  FORMAT  (*  .  FO  .'BAS'D.  @B/_F  .") 

READ  (2.15)  FO,  BAND,  DELTAF 

15  FORMAT  (F8.0,X,F4.0,X,F4.0) 

I  =  E0F(2) 

IF  (FO  .EQ.  0)  FO  =  FOO 

IF  (BAND  .EQ.  0.)  BAND  =  BANDO 

IF  (DELTAF  .EQ.  0)  DELTAF  =  DELTAO 

IF  (FO  »  BAND  *  DELTAF  .EQ.  0)  GO  TO  10 

BANDO  =  BAND 

DELTAO  =  DELTAF 

FOO  =  FO 

NFREQ  =  1 

IF  (DELTAF  .NE.  0.)  NFREQ  =  0.5  +  BAND*FO/( 10Q*DELTAF)+1 

IF  (NFREQ  .EQ.  2»(NFREQ/2))  NFREQ  =  NFREQ  +  1 

RETURN 

END 
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SUBROUTINE  VPLTPT  (12,  F) 


PURPOSE:  This  subroutine  plots  the  patterns  of  a  rectangular 

microstrip  antenna.  j 

PARAMETERS:  12  is  the  current  calculated  to  flow  through  port  2  1 

when  a  load  impedance  (found  in  VLOAD)  is  attached  to  ! 
it.  F  is  the  frequency. 

OPTIONS:  Both  polarizations  in  two  planes,  X-Z  and  Y-Z  are  plotted' 
if  the  "linear  option"  is  chosen.  If  the  "CP"  1 

(circular  polarization)  option  is  chosen,  then  the 
response  of  a  rotating  dipole  is  simulated  for  the  two  ‘ 
aforementioned  planes.  All  scales  are  linear  (not  dB) .  f 
If  the  "individual  normalization"  option  is  chosen,  then 
each  pattern  is  normalized  by  its  own  maximum  over  the 
scan.  The  relative  "cross-pol"  cannot  be  seen  with  this 
option.  At  the  end  of  each  plot,  a  character  must  be 
input.  The  characters  that  are  allowed  have  the 
following  meanings: 

CHARACTER  STANDS  FOR  EFFECT 

S  status  Next  plot  is  same  type  as 

previous  plot. 

M  manual  Ask  for  options  for  next  plot. 

A  automatic  Sam®  -s  "S"  except  no  further 

input  of  options  are  possible 
and  the  plot  is  automatically 
copi ed . 

R  re-plot  Re-plot  last  graph. 

T  terminate  Terminate  all  plotting. 


COMPLEX  ET,  ET1 ,  EP,  EP1,  12  , 

INTEGER  P,  ANS,  NORM,  CP,  FREQ(2) 

COMMON  /DELTA/  DELTA 
COMMON  /OPT/  ANS 
COMMON  /JID/  ID 

COMMON  /II/  AZ,  B,  T.  D,  DIEL.  LOSS,  SIGM,  P.  XP,  YP,  XPP,  YPP,  L, 

*  FO,  NfREQ .  KO,  TWOPI,  PD02A,  PXPOA,  BOA,  BAND,  PI, 

*  ETAO.  K,  GAIN,  DELTAF 

^DIMENSION  EXZX(7|1h  EXZYU21),  EYZX(721),  EYZY(721 ) ,  ID(3), 

DATA  NORM  /0/, INITIAL  /0/ 

IF  (INITIAL  .NE.  0)  GO  TO  400 
IXZ  =  "X-Z;" 

IYZ  =  "Y-Z;" 

10  =  "0;" 

IDASH  =  "-;" 

ISLASH  =  "/:" 

PIBY180  =  PI /I 80. 

PIBY2  =  PI/2. 

PIBYo  =  PI  /  8 
DO  300  J  ■  1.  8 
ANGLE  =  (J-1)  •  PIBY8 
C(J)  =  COS (ANGLE) 

300  S(J)  =  SIN (ANGLE; 

400  IF  (NORM  .NE.  0)  GO  TO  77 

PRINT  *  I 

PRINT  •;  »‘‘riTl*t',lTlT!*m,,!*l*l,,l''l*lT!*!*l*l"  1 
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PRINT  », 

PRINT  »,  "CP  OR  LINEAR  (TYPE  C  OR  L)  ", 

READ  75,  CP 

IF  (CP  .EQ.  "C")  GO  TO  76 

FACTOR  =  PIBY1 80 

FACT  =  1.0 

DENSE  =  91 

LIMIT  =  181 

PRINT  », 

PRINT  •,  "INDIVIDUAL  NORMALIZATION  (Y  OR  N)  ", 

READ  75,  IQ 

75  FORMAT  (A1) 

GO  TO  77 

76  FACTOR  =  PIBY180/4. 

FACT  =0.25 

DENSE  =  361 
LIMIT  =  721 

77  DO  3  J  =  1,  LIMIT 

IF  (CP  .EQ.  "C")  IPSI  =  MOD( J-1 ,8)  +  1 
ANGLE  =  (DENSE-J)  *  FACTOR 

,0.) ,  XP,  YP) 


N/S 


N/S 

N/S 


CALL  VPAT  (ANGLE,  0..  ET,  EP,  DELTA,  (1, 

IF  (P  .EQ.  1)  GO  TO  1 
CALL  VPAT  (ANGLE,  0.,  ET1 ,  EP1 ,  DELTA,  12,  XPP,  YPP) 

EP  =  EP  +  EP1 
ET  s  ET  +  ET  1 

1  IF  (CP  .EQ.  "C")  GO  TO  111 
EXZX(J)  =  CABS(ET) 

EXZY(J)  =  CABS(EP) 

GO  TO  112 

111  EXZX(J)  =  CABS(EP#C(IPSI )+ET*S(IPSI ) ) 

112  CALL  VPAT  (ANGLE,  PIBY2 ,  ET,  EP,  DELTA,  (1.,0.),  XP,  YP) 
IF  (P  .EQ.  1)  GO  TO  2 

CALL  VPAT  (ANGLE,  PIBY2 ,  ET1 ,  EP1 ,  DELTA,  12,  XPP,  YPP) 
EP  =  EP  +  EP1 
ET  =  ET  +  ET1 

2  IF  (CP  .EQ.  "C")  GO  TO  22 
EYZX(J)  =  CABS(EP) 

EYZY(J)  =  CABS(ET) 

GO  TO  3 

22  EYZY(J)  =  CABS(EP*C(IPSI )  +  ET*S(IPSI)) 

IPSI  =  MOD (J-1, 8)  +  1 

3  CONTINUE 
CALL  URESET 
A1  =  0. 

A2  =  0. 

A?  =  0. 

A4  =  0. 

DO  5  J  =  1,  LIMIT 

A1  =  AMAXUA1,  EXZX(J)) 

IF  (CP  .NE.  "C"T  A2  =  AMAX1 (A2,  EXZY(J)) 

IF  (CP  .NE,  "C")  A3  =  AMAX1 (A3,  EYZX(J)) 

5  A4  =  AMAX1 (A4 ,  EYZY(J)) 

A  =  AMAX1  (A1 ,  A2,  A3,  A4) 

CALL  UPEN  (0.,  0.) 

7  CALL  UERASE 


CALL  UDAREA  (0.4,  5.119,  0.4,  5.119) 
CALL  UWINDO  (-1.,  1.,  -1.,  1.) 

8  CALL  UCRCL$  ?0.}  0.,  0.2  •  J) 


vnuu  uyiivuu  \  v  •  •  v  • 

CALL  USET  ("POLAR") 

DO  9  J  =  1.  36 
ANGLE  =  J  *  10 
CALL  UMOVE  (0.2.  ANGLE) 
9  CALL  UPEN  (1..  ANGLE) 

DO  10  J  =  1,  4 
ANGLE  =  J  *  90 
CALL  UMOVE  (0.,  0.)  % 

10  CALL  UPEN  (0.2.  ANGLE) 
CALL  USET  ("LINE") 
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AA  =  A 

IF  (IQ  .EQ.  "Y")  AA  =  A1  % 
CALL  UMOVE  (EXZX( 1 )/AA,90 . ) 

DO  11  J  s  1,  LIMIT 

ANGLE  =  (DENSE  -  J)  «  FACT  v 

11  CALL  UPEN  (EXZX(J)/AA,  ANGLE) 
IF  (IQ  .EQ.  "Y")  AA  =  A3  „ 
CALL  UMOVE  (EYZY( 1 )/AA,90. ) 

DO  12  J  =  1,  LIMIT 

ANGLE  =  (J  +  DENSE  -  2)  «  FACT 

12  CALL  UPEN  (EYZY(J)/AA,  ANGLE) 
IF  (IQ  .EQ.  "Y")  AA  =  A2 

IF  (CP  .EQ.  "C")  GO  TO  120 
CALL  USET  ("RECT") 

CALL  UMOVE  (0.8,  -0.8) 

CALL  UPEN  (1.0,  -0.8) 

CALL  USET  ("POLAR'’) 

CALL  USET  ("DASHLINE") 

CALL  UPSET  ("SETDASH".  5212.) 
CALL  UMOVE  (EXZY( 1 )/AA,90. ) 

DO  13  J*  1,  181 
ANGLE  =  (91  -  J) 

13  CALL  UPEN  (EXZY(J)/AA,  ANGLE) 
IF  (IQ  .EQ.  "Y")  AA  =  A4 
CALL  UMOVE  (EYZX(1 )/AA,90. ) 

DO  14  J  =  1,  181 

ANGLE  =  (J  +  89) 

14  CALL  UPEN  (EYZX(J)/AA,  ANGLE) 
CALL  USET  ("RECT") 

CALL  UMOVE  (0.8,  -0.9) 

CALL  UPEN  (1.0,  -0.9) 


CALL  UPEN  (1 
CALL  UDAREA  < 
CALL  UPRINT  l 
CALL  UPRINT  l 
CALL  UPRINT  i 
CALL  UPRINT  i 
GO  TO  121 


(0(,  6.,  0..  5.11?) 
0.?6,  -0.6$8,  10) 
9-76,  -0.768,  10 


0.76,  -0.688,  IDASH) 
(0.76,  -0.768,  ISLASH) 


120  CALL  USET  ("TEXT") 

CALL  USET  ("RECT") 

CALL  UDAREA  (0.,  6..  0.,  5.119) 

121  CALL  UPRINT  (-1.0.  0.,  IYZ) 

CALL  UPRINT  (+0.76,  0.,  IXZ) 
CALL  UPRINT  (-0.40,  -1.0,  ID) 


CALL  USET  ("LINE") 
ENCODE  (11, 16. FREQ 


16  FORMAT 


(11 ,16. FREQ) 
("F  =  ,F6. 1 , 


CALL  UPRINT  (0 .60 ,-1 .6 ,FREQ) 
IF  (NORM  .NE.  2)  GO  TO  15 
CALL  UFLUSH 
PRINT  *  "  \  >  & 

CALL  UPAUSE 
RETURN 


,  UREAD 

(-1., 

,  -1., 

KQ. 

KQ 

•  EQ. 

"S" 

NORM 

=  1 

KQ 

.EQ. 

"M" 

NORM 

=  0 

KQ 

.EQ. 

"A" 

NORM 

=  2 

,KQ 

.EQ. 

"R" 

GO  TO 

i  7 

KQ 

.EQ. 

"T" ' 

AN  S  = 

3 

1.,  FLAG) 


CALL  UERASE 

RETURN 

END 
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SUBROUTINE  VPLTS  (F,  Z11,  Z12,  Z22,  IACC) 


PURPOSE:  This  routine  plots  the  "s"  parameters  of  a  two  port 
microstrip  antenna  using  the  input  values  of  the 
open  circuit  parameters.  It  uses  the  software  of  the  GCS 
system. 

PARAMETERS:  F  is  the  frequency  in  MHz. 

Z11,  Z12,  and  Z22  are  the  respective  open  circuit 
parameters  of  the  microstrip  antenna  (input). 

IACC  is  set  equal  to  zero  when  data  is  being 

acculmulated  by  the  subroutine  for  later  plotting. 
It  is  set  to  unity  when  the  acculmulated  data  is 
to  acctually  be  plotted. 

OPTIONS:  At  the  end  of  each  plot,  one  may  input  either  a  blank 
character  or  an  "R".  Tne  latter  causes  the  system  to 
replot  the  graph. 


COMPLEX  S1 1(100),  S12 ( 100) ,  S22(100),  Z11,  Z12,  Z22 
DIMENSION  ID(3) 

INTEGER  ONEONE.  ONETWO,  TW0TW0,  S 
COMMON  /JID/  ID 

COMMON  /II/  A.  B,  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP ,  YP ,  XPP ,  YPP,  L, 

*  F0,  NFREQ,  K0,  TWOPI,  PD02A,  TPXPOA,  BOA,  BAND,  PI, 

*  ETAO,  K.  GAIN,  DELTAF 
COMMON  /DELTA/  DELTA 

S  =  "S: " 

IX  =  "X;" 

10  =  "0;" 

IS  =  "+;" 

IC  s  "*;" 

ONEONE  =  "  11;" 

ONETWO  =  "  12;" 

TWOTWO  =  "  22;" 

ID (3)  =  ";" 

IF  (IACC  ,E0.  1)  GO  TO  5 

CALL  VZTOS  (Z11,  Z12.  Z22,  S11(K),  S12(K),  S22(K)) 

IF  (K  .GT.  1)  GO  TO  3 
Q  =  1  /  DELTA 
CALL  DATE  (ID(1)) 

CALL  TIME  (ID (2)) 

WRITE  (3,1)  A,  B,  T.  D,  XP.  YP.  XPP,  YPP,  DIEL,  LOSS, 

*  SIGM,  FO,  DELTAF,  Q,  GAIN,  ID(1),  ID(2) 

1  FORMAT("1 ",6X , "T  WO  PORT  RECTANGULAR  MICRO 

1  S  T  R  I  P" ,/ ,7X ," -  -  -  -  -  -  - 

2-  -  -  -  -  -  -  ,/,"0",6X ."Dimensions  (a  X  b  X  t)  . 

3  .  (*,  F4.1,"  X  ",F4.1,"  X  ", F4.2  ")  cm",/,7X,"F 

4Eed  data:  (1)  Width  .  ",F4.2,"  cm 

5",/,l8X,"(2)  Port  one  .  (",F4.1,",",F4.1 ," 

6)  cm",/,18X,"(3)  Port  two  .  (",F4.1,",", 

7F4.1,")  can", /,7X, "Electrical  parameters:  (1) 

oDi electric  .  ",F5.2,/,30X,"(2)  Loss  tangent 

§  ...  ",F7.5,/,30X,"(3)  Conductivity  ...  ",F5.1,"  KMho 
•s/Cm" ,/ ,7X , "Center  frequency  . 

1  .  ",F7.1,"  MHz", /,7X, "Frequency  increment  .. 

2  .  "  ,F7.1,"  MHz", /,7X, "Quality  fact 

3or  .  ",F5.1./,7X,"Gain  . 

4  .  " ,F4. 1 ,"  dB" , / ,7X , "Date/time  . 

§ . . . . .  ",2A10,/."0".6X,"Frequenc' 

oY" , 1 IX ,"S" , loX ,"S" , loX ,"S" ,/ ,9X ,"(MHz) ", 14X,"1 1" , 15X,"12" , 15X,"12 
7" ,/ ,"0") 

3  WRITE  (3,4)  F,  S11(K).  S12(K),  S22(K) 

4  FORMAT  (fX,F7.1,6X,3(*(",F6.3,","F6.3,")",2X)) 

RETURN 
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CALL  URESET 
CALL  UERASE 
N02  =  NFREQ/2  +  1 
N021  =  N02  +  1 
CALL  USET  ("LINE") 


CALL  UWINDO  7-1 .. 1..-1. ,1.) 

CALL  UCRCLE  i  0 ! '  ’  o!  I  ’  1?  5  5'"9> 

CALL  UMOVE  (0.  A.)  ’ 


CALL  UAHC*'(i:;K,90.) 

CALL  UARC  ( 1.1-1 1,90.) 

CALL  UCRCLE  (6. 5, 6. 0.0. 5) 
CALL  UMOVE  (-1.oIo.oJ 
CALL  UPEN  (1. 0,0.0) 

CALL  USET( "DASHLINE") 

CALL  UPSET ("SETD ASH", 54 14.) 
CALL  UCRCLE  (0.,0.,0,5) 

CALL  USET("LINE*) 

Aif  T  ITnnm  / 


CALL  USET  ("NX")' 

CALL  USET  ("SOFT") 

CALL  UPSET  ("HORIZONTAL".  0.03) 

CALL  UPSET  ("VERTICAL",  6.03) 

DO  6  J  =  1,  NFREQ 

IF  (J  .EQ.  N02)  CALL  USET  ("N*") 

IF  (J  .EQ.  N021 )  CALL  USET  ("NX") 

6  as  *f^T{i:ifo«feiw(sii(j))) 

CALL  USET  ("NO") 

DO  7  J  =  1,  NFREQ 

W  *EQ.  N02)  CALL  USET  ("N«") 

_  IF  (J  .EQ.  N021 )  CALL  USET  ("NO") 

CALL  UPSET  ("VERTICAL",  0.06) 

CALL  UPSET  ("HORIZONTAL",  0.04) 

DO  8  J  =  1,  NFREQ 

H  *IQ-  N02 )  CALL  USET  ("N»") 

0  IF  (J  .EQ.  N021 )  CALL  USET  ("N+") 

CALL  USET  ("HARD")  ' 

as  upriSt  {8:11; 

as  SSI  8:8:  3:83:  i 

CALL  UPRINT  (0.85,  -0.90,  ONETWO) 

CACL  UPRINT  (0.85,  -1.,  fwOTWO) 

CALL  UDARGA  ( 0  • •  5  •  5  •  0  • *5*5) 

CALL  UPRINT  (-0  47,-1 .6, ID) 

CALL  UREAD  (-1..  M  ’  l<5,  i\,  F) 

IF  (IQ  ,NE.  "  "J  GO  fo  5  * 

CALL  UERASE 

RETURN 

END 


DO  7  J  =  1, 
IF  (J  .EQ. 
IF  (J  .EQ. 


CALL  UPRINT 
CALL  UPRINT 
CALL  UPRINT 
CALL  UPRINT 
CALL  UPRINT 
CALL  UPRINT 


nooooonnooooono 


SUBROUTINE  VPLTZ  (Z,  F,  IACC) 


PURPOSE:  This  program  accumulates  impedance  data,  outputs  in  tab-  1 
ulated  form,  and  then  plots  in  one  of  three  optional  ways:  1 
(a)  Smith  chart,  (b)  Rectangular  G-B  plot,  or  (3)  magnitude1 
of  reflection  coefficient  vs  frequency. 

PARAMETERS:  Z  is  the  complex  impedance  to  be  plotted. 

F  is  the  frequency  at  which  Z  was  determined. 

IACC  is  2,  data  is  accumulated  for  future  plotting,  ‘ 
while  if  IACC  is  1,  the  data  is  plotted.  1 


COMPLEX  Y ,  GAMMA,  Z 
INTEGER  P,  STRING(4) 

DIMENSION  ID(3) ,FREQ( 100 ) ,G( 100 ) ,B( 100 ) ,IF0 (2) ,IINC(2) ,IDIM(2) 
COWON  /JID/  ID 
COMMON  /LDID/  STRING 

COMMON  /II/  A.BB.  T,  D,  DIEL,  LOSS,  SIGM,  P,  XP,  YP,  XPP,  YPP,  L, 

•  FO,  NFREQ,  KO,  TW0PI,  PD02A ,  TPXPOA,  BOA,  BAND,  PI, 

•  ETAO,  K,  GAIN,  DELTA* 

COMMON  /DELTA/  DELTA 

ID(3)  =  ";" 

IF  (IACC  .Eq.  1)  GO  TO  10 
N02  =  NFREQ/2  +  1 
N021  =  N02  +  1 
IF  (K  .GT.  1)  GO  TO  3 
Q  =  1  /  DELTA 
CALL  DATE  pDjljj 


CALL  TIME  _ 

ENCODE  (17,  202,  IDIM)  A,  BB 
202  FORMAT  (F5.2  "  X  ",F5.2  *  CM;") 

ENCODE  (17.  200,  IFO)  FO 

200  FORMAT  ("FO  =  *,F6.1,"  MHZ;") 

ENCODE  (16.  201,  IINC)  DELTAF 

201  FORMAT  ("INC  =  *  F5.1,"  MHZ;") 

WRITE  (3,D  A,  BB,  T,  D.  XP,  YP 

IF  (P  .EQ.  2)  WRITE  (3,999 )  XPP.  YPP,  STRING 
WRITE  (3,1000)  DIEL, LOSS, SIGM, Fu , DELTAF, Q, GAIN , ID (1 ) ,ID(2) 
1  FORMAT(*1",6X,"T  WO  PORT  RECTANGULAR  M 


N/S 

N/S 

N/S 


i  a-  vmiA  i  \  i  •  ua  «  i  n 

1  S  T  R  I  P"  ,/^7X," 


I  C  R  0 


2 —  -  -  -  -  -",/,"0",6X, "Dimensions  (a  X  b  X  t)  .... 

3 . .  (*,  F5.2,"  X  ",F5.2,"  X  ".F4.2  ") 

4Eed  data:  (l)  Width  .  *,F4.2," 


cm" ,/  ,7X  ,"F' 


cm 


1000 


",/,18X,"(2)  Port  one  .  (",F5.2,",",F5.2," 

o)  cm") 

999  F0RMATO8X  "(3)  Port  two  .  (", 

1F5.2  ",",F5.2,")  Cm",/, 

2JX^oad  description:  .  ", 

FORMAT (  7X, "Electrical  parameters:  (1) 

SDielectric  .  ",F5.2,/,30X,"(2)  Loss  tangent 

9  ...  " ,F7.5,/ ,30X ,"(3)  Conductivity  ...  ",F5.1,"  KMho 
*s/ Cm", /,7X, "Center  frequency  . 

1  . ",F7.1,"  MHz", /,7X, "Frequency  incr«nent  .. 

2  .  "  ,F7 . 1 ,"  MHz", /,7X, "Quality  fact 

3or  . !.  ".F5.1,/,7X,"Gain  . 

4  .  ",F4.1,"  dB",/,7X, "Date/time  . 

5  . . .  ",2A10,/,"0",6X,"Frequenc' 

6Y".11X,"Z",/,9X,"(MHz)", 

7/,^0") 

WRITE  (3,4)  F,  Z 

FORMAT  (7X,F7.1,6X,"(",F6.2,","F7.2,")") 

FREQ(K)  =  F  -  FO 
Y  =  50. /Z 
G(K)  =  REAL(Y) 


B(K)  =  AIMAG(Y) 
RETURN 

10  CALL  URESET 
PRINT  •  "  • 
PRINT  • 

PRINT  * 

PRINT  * 

PRINT  • 

1 

READ  2, 

2  FORMAT 


„.rr  rrrrrrrrrr  rrrrrrrrrr!„ 

"RECTANGULAR,  SMITH  CHART,  OR  VS  FREQUENCY  PLOT  (R,  S,  ", 
"OR  V)", 


14 


12 


13 


,  KODE 
CAD 


GO  TO 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


(E0F(2) 


(KODE 

(KODE 

(KODE 


W  f  »« 

UERASE 
USET  ("LINE") 
UWINDO  ' 
UDAREA 
UCRCLE  . 
UMOVE  (0 


RETURN 
KKODE  = 
KKODE  = 
KKODE  = 
KKODE 


-1.,1.,-1.,1.) 

0.4,  5.119.  0.4,  5.119) 

0  0.,  1.5 
unu  v  Ei  \  u  •  •  I  •  ) 

UARC  ( 1 . , 1 . , 90 . ) 

UARC  (1..-1..90.) 

CALL  UCRCLE  (0 .5 ,0 .0. 0 .5) 

CALL  UMOVE  (-1 .0,0.0) 

CALL  UPEN  (1. 0.0.0) 

CALL  USET ("DASHLINE") 

CALL  UPSET ( "SETDASH" ,5414.) 

CALL  UCRCLE  (0..0..0.5) 

CALL  USET ("LINE*) 

CALL  USET  ("NX") 

CALL  USET  ("SOFT") 

CALL  UPSET  ("HORIZONTAL".  0.02) 

CALL  UPSET  ("VERTICAL",  0.02) 

DO  6  J  =  1,  NFREQ 
Y  =  CMPLX  (G(J).B(J)) 

GAMMA  =  (1-Y)/(1+Y) 

IF  (J  .EQ.  N02)  CALL  USET  ("N«") 

IF  (J  .EQ.  N021 )  CALL  USET  ("NX") 

CALL  UPEN  (REAL (GAMMA) ,  AIMAG (GAMMA ) ) 
USET  ("HARD") 

UDAREA  (0., 5. 5,0.. 5. 5) 

UPRINT  (-0.40, -1.0, ID) 

USET  ("SOFT") 

USET  ("ITALLICS") 

UPSET  ("HORIZONTAL",  0.03) 

UPSET  ("VERTICAL",  0.05) 

UPRINT  (0.4,  -0.85,  IFO) 

UPRINT  (0.4,  -0.92,  IINC) 

UPRINT  (0.4,  -0.99,  IDIM) 

UERASE 

T  UREAD  (0..  0.,  KODE,  1.,  FLAG) 
(KODE  .EQ.  "R")  GO  TO  10 
TO  9 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
IF  ' 

GO 

RMAX  =  0. 

XABSMX  =  0. 

DO  13  J  =  1,  NFREQ 
Z  =  1  /  CMPLX  (G(J).  B(J)) 

XABSMX  =  AMAX1  (ABS (AIMAG ( Z) ) ,  XABSMX) 

RMAX  =  AMAXKRMAX,  REAL(Z)) 

ABSMAX  =  AMAX1  (RMAX,  XABSMX) 

DEC  =  AL0G10( ABSMAX) 

IDEC  =  DEC 

IF  (DEC  .LT.  0.)  IDEC  =  IDEC  -  1 
DECADE  =  10.  «  IDEC 
IDIGIT  =  ABSMAX  /  DECADE  +  1 
WINDOW  =  IDIGIT  *  DECADE 

CALL  UWINDO  (0.,  2.* WINDOW,  -WINDOW,  WINDOW) 
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-j; «»» •  xtwi’-v.  ■ .  j?- --  -'•* 


TICK  =  DECADE 
IF  (IDIGIT  .LT.  7)  TICK  = 
IF  (IDIGIT  ,LT.  3)  TICK  = 
CALL  UPSET  ("TICX" ,  TICK) 
CALL  UPSET  ("TICY",  TICK) 
CALL  UERASE 

CALL  UDAREA  (0.4,  5.,  0.4 
CALL  UPSET  ("YLABEL","X;" 


7)  TICK  =  0.5  *  DECADE 
3)  TICK  =  0.2  «  DECADE 


vnuu  uunnun  ^ J • t 

CALL  UPSET  ("YLABELV'X;”) 

CALL  UPSET  ("XLABEL" ,"R; ") 

CALL  USET  ("GRID AXIS") 

CALL  USET  ("XBOTH") 

CALL  USET  ("YBOTH") 

CALL  USET  ("POINT") 

CALL  USET ( "OWNSCALE" ) 

CALL  UAXIS  (0.,  2. "WINDOW,  -WINDOW,  WINDOW) 
CALL  USET ("LINE") 

CALL  USET  ("NX") 

CALL  USET  ("SOFT") 

CALL  UPSET  ("HORIZONTAL".  0.02  *  WINDOW) 
CALL  UPSET  ("VERTICAL",  0.02  •  WINDOW) 

DO  130  J  s  1,  NFREQ 
Z  =  1  /  CMPLX  (G(J),  B( J ) ) 

130  CALL  UPEN  (REAL(Z) ,  AIMAG(Z) ) 

CALL  USET  ("HARD") 

CALL  UDAREA  (0.,  5.,  0..  5.) 


!v>.,  J.  ,  u.  .  J.  / 

-1.,  1.,  -1.,  1.) 

-0.47, -1.0, ID) 

CALL  UPAUSE 
CALL  UERASE 
CALL  USET  ("LINE") 

GO  TO  14 
15  CALL  UERASE 

CALL  UDAREA  (0.4.  5.4.  0.4,  5.4) 

CALL  USET  ("GRID AXES") 

CALL  USET  ("XBOTH") 

CALL  USET  ("YBOTH") 

CALL  UPSET ("XLABEL", "FREQUENCY  -  FO:") 

CALL  UPSET  ( "YLABEL" , "MAGNITUDE  OF  GAMMA;") 

CALL  USET  ("AUTO") 

CALL  UAXIS  (FREQ(1),  FREQ (NFREQ ) ,  0.00,  1.00) 

CALL  USET ("SOFT") 

CALL  UPSET("HORIZONTAL" ,0.02  *  (FREQ(NFREQ)  -  FREQ ( 1 ) ) ) 
CALL  UPSET ( "VERTICAL" ,0.02) 

CALL  USET ("NX") 

DO  156  J  =  1.  NFREQ 
Y  =  CMPLX  (GU).B(J)) 

GAMMA  =  (1-Y) /( 1+Y) 

156  CALL  UPEN  (FREQ(J) ,  CABS ( GAMMA ) ) 

CALL  USET  ("SOFT") 


CALL  USET  ("ITALLICS") 

CALL  UDAREA  (0.,  7.49.  0..  5.71) 
CALL  UWINDO  (o.,  1..  0..  1.) 

CALL  UPSET  ("HORIZONTAL*.  0.015) 
CALL  UPSET  ("VERTICAL",  0.03) 


CALL  UPRINT 
CALL  UPRINT 
CALL  UPRINT 
CALL  UPRINT 
CALL  UPAUSE 
CALL  UERASE 
GO  TO  14 
9  RETURN 
END 


0.7,  0. 
0.7,  0. 


ID)  % 

,  IFO). 


v  •  |  v  / 

0.7,  0.44,  IINC) 
0.7,  0.41,  IDIM) 


CHAPTER  3:  EXAMPLES 


This  chapter  contains  examples  of  the  use  of  this  program  to  analyze  a  microstrip  antenna. 
The  specific  case  chosen  was  that  of  a  nearly  square  microstrip  antenna.  This  was  chosen  to 
illustrate  that  by  making  one  side  of  the  antenna  a  small  amount  larger  than  the  other,  and  by 
properly  feeding  and  loading  the  the  antenna  with  a  variable  capacitor  (such  as  a  varactor),  the 
antenna  can  be  switched  from  left  hand  circular  polarization  to  right  hand  circular  polarization. 
This  is  consistent  with  both  theory  and  experiments  carried  out  at  the  University  of  Illinois  and 
presented  in  references  [4]  and  [6].  The  output  listed  below  is  a  copy  of  the  actual  graphical 
data  displayed  on  a  graphics  terminal  by  the  program.  Explanatory  remarks  have  been  added  to 
aid  the  reader. 
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